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Abstract 



O 

o 
o 
o 

We show that at least the quasi-exactly solvable eigenvalues of the 
' Schrodinger equation with the potential V{x) = —{^ cosh. 2x — iM) as well as 

the periodic potential V{x) = (("cos 20 — iM)'^ are real for the PT-invariant 
non-Hermitian potentials in case the parameter M is any odd integer. We fur- 
ther show that the norm as well as the weight functions for the corresponding 
weak orthogonal polynomials are also real. 
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In non-relativistic quantum mechanics, one usually chooses a real (Hermitian) potential 
so as to ensure the real energy eigenvalues of the corresponding Schrodinger equation. About 
three years ago, Bender and others have studied several complex potentials which are 
invariant under the combined symmetry PT and showed that even in all these cases, the 
energy eigenvalues of the Schrodinger equation are real. This seems to suggest that instead 
of Hermiticity, it may be enough to have PT-invariant Hamiltonian so as to have real energy 
eigenvalues. It must however be noted that so far, this is merely a conjecture which is being 
supported by several examples. Besides, even if this conjecture is true, there are several basic 
questions which will have to be addressed before one can take the PT-invariant potentials 
more seriously. 

In the absence of any rigorous proof, it is worthwhile studying new PT-invariant non- 
Hermitian potentials and see if in each and every case the eigenvalues are real or not. 
The purpose of this note is to study this question through two concrete examples. In 
particular, we consider the system described by the non- Hermitian PT-invariant Hamiltonian 

{h = 2m = 1) 

H =p^ ~{Ccos\i2x-iMf , (0.1) 

where the parameter C, is real and parameter M has only integer values. We show that 
at least the quasi-exactly solvable (QES) eigen values @] of this Hamiltonian are real and 
bounded from below. We also show that the corresponding set of polynomials form a family 
of weak orthogonal polynomials and further, the norms for these polynomials as well as 
the corresponding weight functions are real. Using the duality transformation, we further 
show that the QES eigenvalues of the complex, non-hermitian, but PT-invariant periodic 
potential 

y(^) = (Ccos2^-iM)2, (0.2) 

are also real. 

Let us first show that the Hamiltonian ( |0.1| ) is PT-invariant. To that end, let us first 
notice that under T reflection one replaces i by —i while under P reflection one replaces 
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X hj a — X where a/2 is the origin about which one is performing the parity reflection. It 
is then easily checked that the Hamiltonian (p.l| ) is indeed PT symmetric under the parity 
reflection x i7T/2 — x. 
We substitute 

^(x) = e^^^^^^Xx), (0.3) 

in the Schroodinger equation Htp = Eip with H as given by Eq. ( p.l[ ) and obtain 

(j)"{x) + 2iC sinh 2x(j)'{x) + [{E - + C^) - 2i{M - 1)C cosh 2x] (pix) = . (0.4) 

Ordinarily, the boundary conditions that give quantized energy levels are tplx) ^ as \ x \^ 
oo on the real axis. However, in the present case, we have to continue the eigenvalue problem 
into the complex-x plane ||^. On putting, x = u + iv where u,v are real, it is easy to see 
that for u > the above boundary condition is satisfied so long as — tt < v < —7t/2 [mod it) 
while for m < it is satisfied if — 7r/2 < i; < {mod tt). 
On further substituting 
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z = cosh2x - 1; = £ ^ ) ' , (0.5) 

n=0 V ^ / 

we obtain the three-term recursion relation (n > 0) 

Rn+2{E) - [n^ + 4(s - iOn + As^ + E - + C + 2i{M - !)(] R^{E) 

- AC[M + l-2s-n]n{n-l)Rn-2{E)=Q, (0.6) 

provided 2s^ = s i.e. either s = or s = |. Thus we have two sets of independent 
solutions; one for s = and other for s = |. Note that the parameter s is not contained 
in the potential and this is perhaps related to the fact that for any odd M (> 1), the QES 
solutions corresponding to both even and odd number of nodes are obtained. 

^From Eq. ( p.6| ) we observe that the even and odd polynomials Rn{E) do not mix with 
each other and hence we have two separate three-term recursion relations depending on 
whether n is odd or even. In particular, it is easily shown that the three-term recursion 
relations corresponding to the even and odd n cases respectively are given by (n > 1) 



Pn{E) - [4n2 + 8n{s - < - 1) + 4^^ - 8s + 4 + 6iC + - (M - iCY\ Pn-i{E) 

+ 8<(n - 1) {2n - 3) (M + 3 - 2s - 2n) P„_2 = , (0.7) 

Qn{,E) - [Au^ + 4n(2s - 2< - 1) + 4s2 - 4s + 1 + 2< + E - (M - iCf] Qn-i{E) 

+ 8<(n - l)(2n - 1) (M + 2 - 2s - 2n) g„_2 = , (0.8) 

with Po{E) = l,QoiE) = 1. 

We find that for odd integer values of M, there are M quasi-exactly solvable eigen values 
and all of them are real provided |C| < Co where Cc is a function of M. At C = Cc two highest 
eigen values become degenerate, while for C > Cc some of the eigen values are complex. In 
particular, if M = 2/c + 1 {k = 0,1, ...), then k + 1 and k levels with real energy eigenvalues 
are obtained from the zeros of the orthogonal polynomials Pk+i{E) and Qk{E) respectively 
which we term as the critical polynomials. Further, all higher P{Q) polynomials exhibit 
factorization property i.e. 

Pfc+„,+i(E) = Pfc+i(E)P„(E) 
Qk+niE) = QkiE)Q^{E) , (0.9) 

where Po{E) = 1 = Qq{E). It is interesting to note that even though the polynomials 
satisfying recursion relations ( p.7|) and (|0.8| ) are complex, all critical polynomials (i.e. Pk+i 



and Qk) are real in case M = 2k + 1. 

We now consider several odd integer values of M and show explicitly that the QES eigen 
values are real for ( < (c- Let us start with M = 1, in which case Pi{E) = E — 1 + ("^ 
is the critical polynomial and the only QES level is E = 1 — ('^ which is always real. It is 
not very clear if this is the ground state energy of the system or not even though for the 
corresponding real potential, for M=l one indeed obtains the ground state energy 0. 

For M = 3 the critical polynomials are, 

P2 = £2 + 4£ + 16C^ 

Qi=S + 4, (0.10) 
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where 



E = E-M'^ + C. 



(0.11) 



The energy levels are 



E 



5-e 



E 



7 - - 2^1 - 4(2 



E 



7 - + 2^/1 - 4(2 



(0.12) 



of which the first energy is obtained from Qi = while the other two are obtained from 
P2 = 0. For < (^2 _ ]^^^ these levels are real and at C = Cc highest two levels are 
degenerate. 

Proceeding in the same way, we have considered the cases of M = 5, 7, 9 but unfortunately 
in these cases (except when M=5 and Q2 = 0) one has to solve cubic or higher order equations 
so as to obtain the corresponding energy eigenvalues. We have calculated these eigenvalues 
in case = 0.01, 0.02 and 0.025 (all of which are less than the corresponding and they are 
given in Tables I, II and III. In each case, we also give the expressions for the corresponding 
critical polynomials as well as (c- 

We have studied the properties of the weakly orthogonal polynomials Pn{E) and Qn{E) 
and we find that many of these are similar to the Bender-Dunne polynomials 0. In partic- 
ular, since in both the cases the recursion relations are similar to those given by Eq. ( p. 71) , 
hence for all values of the parameters M and s, they form an orthogonal set. Secondly, 
the wave function iIj{x,E) is the generating function for the polynomials Pn{E) as well as 
Qn{E). Thirdly when M is a positive integer, both of these polynomials exhibit factorization 
property whose precise form depends on whether M is even or odd. The norms and weight 
functions for these polynomials can be calculated using the same procedure as described in 
0] and they turn out to be complex. 

We now show that there is an alternative way |]^ to express this system as an algebraic 
quasi-exactly solvable system and that in this case the norms and the weight functions of 
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the corresponding orthogonal polynomials are real. 

Following Finkel et al. p|, we change the variable, z = e^^ and make the gauge transfor- 
mation 

^{z)=z'-^e'-^^'+-^\ (0.13) 

in the Hamiltonian ( p. I]) to map it to a differential operator, Hg ( called as gauge Hamilto- 
nian). This gauged Hamiltonian, Hg can be expressed in terms of the generators of S1(2,R) 
by 

Hg{z) = -AJl - 2iCJ+ + 2<J_ - c* , (0.14) 

where c* = — + while the generators ^+,-,0 are given by 

d Ti d 

J- = g;^ ^o = .^-2^ J_ = .^^-n.. (0.15) 

Now on substituting 

^0 (2<)"n! 

in the Schrodinger equation, Hgcj) = Ecp, one finds that Rn{E) satisfy a three-term recursion 
relation given by 

Rn+liE) = {E- hn)Rn{E) - anRn-liE), n > , (0.17) 

where 

a„ = -An{M - n)C , 

6„ =4n(M- 1 -n) + 2M- 1 -C^ (0.18) 

Thus we see that (with Rq = 1) the same system ( p.l|) forms another family of weakly 
orthogonal, monic polynomials. The first few polynomials generated by the recursion relation 
( p.l7| ) are given by 
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Ri = E-bo 

R2 = E^- (bo + bi)E + ai + b^bo 

R3 = E^ - (bo + bi + b2)E^ + (60&1 + &0&2 + + oi + aa)^ 

- (60^1^2 + 0162 + a2&o) 
R^ = {E-bs)R3-asR2, (0.19) 

where a„ and 6„ are given by the Eq. ( |0.18| ). 

These polynomials are also of the Bender-Dunne type and exhibit similar properties as 
the Bender-Dunne polynomials. For M = positive integer (say k,) these polynomials exhibit 
factorization property given by, 

Rk+n = RkRn n>0. (0.20) 



Note that at M = A;, = and hence Eq. ( p.l7| ) reduces to two term recursion relation. 



Further, when M is odd (say M=2k+1) then not surprisingly we find that the corresponding 
critical polynomial factorizes as the product of Pk+i and Qk, i.e. 

R2k+iiE) = Pk+^{E)Qj,{E) , (0.21) 

where Pk, Qk are as defined by Eqs. (|0.7]) and (p.8|). 

Here the polynomials Rn correspond to the non-QES part of the spectrum and they 
satisfy the three-term recursion relation 

Rn+l = {E ~ ^M+n)Rn — O-M+nRn-l (0.22) 

Hence these polynomials (with = 1) also form a set weakly orthogonal polynomials. 

The norms (squared) of the orthogonal set of polynomials Rn{E) can be easily determined 
using the recursion relations of these polynomials and one finds 

n 

7n = (-l)"n«- (0-23) 

i=l 



where are given by eq. ( p. 18] ). It is clear that the norms for these polynomials, even 



though real, are not always positive. The norms for the even polynomials are positive and 



on the other hand the norms for the odd polynomials are negative so long as n < M while 
for n > M, 7„ vanishes. Similarly, the weight functions, uj{E) for these polynomials are also 
real, though not positive definite. This is easily seen using the proposition 2.1 in and the 
fact that Ofc are negative for 1 < A; < M. 

Before finishing this note, we show that the QES eigenvalues of a related complex non- 
Hermitian but PT-invariant periodic potential are also real. To that end, it is worth noting 
that recently Krajewska et. al. JlOl have discussed the consequences of the anti-isospectral 
transformation (also termed as duality transformation) in the QES problems. In particular, 
they have shown that under the transformations x — > = y, if a potential V{x) goes to 
V{y) then the QES levels of the two are also related. In particular, they have shown that 
if M levels of the potential V{x) are QES levels with energy eigenvalues and eigenfunctions 
Ej:{k = 0, 1, ■ ■ ■ , M — 1) and ipki^) respectively then the energy eigenfunctions of V{y) are 
given by 

Ek = -Eu-i-k, fpkiy) = ipM-i-kiix) (0.24) 

We shall now apply this anti-isospectral transformation to the system described by Eq. 
( p. ID and obtain another example of a non-Hermitian but PT-invariant system with real QES 
energy eigenvalues. On applying the duality transformation x iO to the Schroodinger 
equation Hip = Eip with H as given by Eq. (|0.1j ) we obtain the following Schrodinger 
equation for the PT-invariant DSG potential 

Using Eq. (p.24| ), the QES levels for this system can be obtained for odd integer values 
of M from the QES levels of the dual system which we have already discussed in details. 
For example, for M = 1, the only QES level is, 

^ = -(l-C^), ^ = 6+4^*^26^ (0.26) 

which is real for any value of C, . On the other hand for M = 3 the three QES levels are 
given by 
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{Ccos2e -iMf 



ip{e) = Eipie), E = -E. (0.25) 



-2< - I V^l - 4(2 - l| cos 2^ 



^4 cos 26* 



E = -7 + + 2i/l - 4(2 , V'= -2< + |i/l-4C2 + l|cos2^ 



Similarly one can easily obtain the QES levels for the cases, M = 5, 7, 9 ■ ■ ■ etc by replacing 
E by —E in the corresponding expressions for the complex DSHG case. 

One can now study the Bender-Dunne polynomials of this dual system and it is easy to 
see that most of the discussion of the previous system goes through in this case. 

When M is even integer then the eigenvalues are complex conjugate pair. Further, 
one can show that whereas PT symmetry remains unbroken when M is odd integer, it is 
spontaneously broken when M is even integer [pJ]j . 

Summarizing, we have seen that at least the QES eigenvalues of the non-Hermitian PT- 
invariant potentials ( |0.1j ) and ( p.25| ) are real and bounded from below. It would be nice if 
one can also compute the non-QES spectrum of these potentials and show that even these 
energy eigenvalues are real. 

Acknowledgment: We are grateful to Prof. Carl M. Bender for pointing out the 
PT-invariance of the potentials discussed in this paper. 
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Table I : QES levels for the case M=5 [Cc ^ 08757] 



Energy Levels 


C = -01 


C = -02 


= .025 


Ep 


9.00331480 


9.00659241 


9.00821724 


Eq 


9.00334818 


9.00672620 


9.00842644 


Eq 


20.97665182 


20.95327381 


20.94157356 


Ep 


21.10018724 


21.20836270 


21.26599253 


Ep 


24.86649797 


24.72504489 


24.65079023 



The critical polynomials in this case are given by, 

P3 = + 20£:^ + £(64 + 64C^) + 768C^ 

^2 = £^ + 20£ + 64 + 16C^ (0.28) 
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Table II :QES levels for the case M=7 ^ 04435] 



Energy Levels 


C = -01 


C = -02 


= .025 


Eq 


13.00199970 


13.00399883 


13.00499708 


Ep 


13.00199971 


13.00399878 


13.00499918 


Ep 


33.01123852 


33.02228660 


33.02773885 


Eq 


33.01140573 


33.02295762 


33.02878902 


Eq 


44.95659457 


44.91304361 


44.89121290 


Ep 


45.21134926 


45.46170007 


45.60896931 


Ep 


48.73541251 


48.43201450 


48.25829367 



The critical polynomials for this case are, 

P^ = S^ + 56£^ + (784 + 160C^)^^ + (2304 + 6528C^)^ 
+ 2304C^ + 55296C^ 

= £3 ^ 56^2 + (784 + 64C^)^ + 2304 + 1536C^ (0.29) 



12 



Table III : QES levels for the case M=9 ^ 02675] 



Energy Levels 


e = .01 


e = -02 


= .025 


Ep 


17.00096869 


17.00182564 


17.00225406 


Eq 


17.00153729 


17.00295705 


17.00366471 


Eq 


45.00607980 


45.01268446 


45.01599526 


Ep 


45.00751350 


45.01552516 


45.01952952 


Ep 


65.02149678 


65.04164257 


65.05150021 


Eq 


65.02298207 


65.04558083 


65.05703252 


Eq 


76.92940084 


76.85877767 


76.8233075 


Ep 


77.37417833 


77.91167617 


78.38610328 


Ep 


80.5458427 


79.92933046 


79.41561294 



The critical polynomials for this case are, 

Ps = + 120£^ + (4368 + 320C^)^^ + (52480 + 30280C^)^^ 

+ (147465 + 806912C^ + 16384C^)^ + 655360C^ + 5898240C^ 
Q4 = + 120^:3 + (4368 + 160C^)^:^ + (52480 + 11648C^)£: 

+ 147456 + 182272C^ (0.30) 
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